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Abstract 

Properties of scalar-isoscalar mesons are analysed in an unitary model 
using separable interactions in three decay channels: tttt , KK and an effec- 
tive 2tt2tt. We obtain different solutions by fitting various data on the tttt 
and KK phase shifts and inelasticities including the CERN-Cracow-Munich 
measurements of the 7r _ pj — ► 7r + 7r~n reaction on a polarized target. Analy- 
tical structure of the meson-meson multichannel amplitudes is studied with a 
special emphasis on the role played by the S'-matrix zeroes. S'-matrix poles, 
located in the complex energy plane not too far from the physical region, are 
interpreted as scalar resonances. We see a wide /o(500) , a narrow /o(980) 
and a relatively narrow /o(1400) . In one of our solutions a resonance at 
about 1700 MeV is also found. Total, elastic and inelastic channel cross sec- 
tions, branching ratios and coupling constants are evaluated and compared 
with available data. We construct an approximation to our model and show 
that the Breit-Wigner approach has a limited phenomenological applicability. 



1 Introduction 

Full classification and identification of scalar mesons is not yet well established 
[[jj . Many theoretical and experimental efforts have been recently made for a better 
understanding of the scalar mesons as it can be seen, for example, in references 
given in @- From QCD one expects presence of some scalar ( J PC = ++ , 1 = 0) 
glueballs which can be mixed with ordinary qq scalar states ||. The lowest scalar 
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glueball masses predicted by lattice QCD calculations are in the range 1500 MeV to 
1700 MeV U|, |5|. There are now lively discussions about the nature of scalar mesons 
/o(1500) and /j(1710) and their possible mixture with scalar glueballs [|], ||. High 
statistics experiments on meson productions such as pp annihilation, nN scattering 
on unpolarized and polarized targets, central production in pp scattering, J/\I/ or 
other heavy meson decays and 77 collisions have been performed. Their analyses 
give some new evidence for the existence of five scalar-isoscalar mesons /o(400— 1200) 
or a, /o(980) , /„(1370), /„(1500) and / (1710) §. 

In we have analysed new solutions for the scalar-isoscalar tttt phase shifts 
[§] together with previous KK results in the framework of a three coupled channel 
model based on an extension of the two-channel model of PJ. Separable potentials 
were used to describe interactions in the tttt , KK and an effective 2tt2tt channels. 
Let us recall that use of separable potentials leads to an analytical solution for the 
S-matrix which helps to check and understand the results. It furthermore represents 
one of the easiest way to handle non-local interactions and so it should be adequate 
to describe, with a minimum number of parameters, strong energy dependence in 
several coupled channel reactions as it seems to be the case here. Simple analytical 
form of separable interactions is then phenomenologically very useful in fitting the 
tttt and KK experimental data from the corresponding thresholds up to 1.8 GeV. 

The third effective 2tt2tt channel was introduced to take into account a strong 
four-pion production observed in different experiments. In these data some evidence 
for four pion clustering into 00 pair, coming from strong interaction between two 
pions was found (see for example [fL0[|). Our 2tt2tt effective channel, called a a , can 
however represent also other possible clusterings such as pp, ai(1260)7r, 7r(1300)7r 
and ujuo. One should not mix the effective threshold mass with twice the mass of 
the /o(500) or a resonance which is seen in the tttt channel. 

The parameters of the model were determined by a fit to two sets of tttt phase- 
shifts and inelasticities obtained in a recent analysis of the CERN- Cracow-Munich 
measurements of the is~p\ — ► TT + TT~n reaction on a polarized target ||. It was 
stressed in || that the a\ exchange gives an important contribution to the ir~p — > 
7T + 7T~n reaction amplitudes. Recently Achasov and Shestakov have also come to 
the conclusion that the a\ exchange plays an important role in the reaction n~p — > 
7T 7r°n on unpolarized target [IT]. This conclusion may be experimentally verified by 
measurements on polarized target. The tttt data |J, covering the 600 to 1600 MeV 
tttt invariant mass range, were completed in the lower energy region by the tttt phase 
shifts of |[12| , [T3|j and [[14]]. Further constraints were imposed by using the KK phase 



shifts from the KK threshold up to 1530 MeV JO 

We found a relatively narrow (90-180 MeV) scalar resonance /o(1400 — 1460). 
Our analysis of previous CERN-Munich unpolarized target data [|16j predicted a 



much broader (r ~ 500 MeV) state. We have also obtained a very wide (r « 500 
MeV) / (500) resonance and the well established narrow / (980) (r w 60-70 MeV). 
Our model allows a theoretical study on the origin of resonances by switching off 
the interchannel couplings. In all the solutions found in the KK interaction was 



2 



repulsive or not attractive enough to create, by itself, a KK bound state. 

In this paper we built two new solutions in which the KK state is bound in the 
uncoupled case. Our solutions are also characterized by presence or absence of a era 
bound state. Gradual increase of interchannel couplings allows one to link the S- 
matrix poles from the uncoupled to fully coupled case. In our approach experimental 
resonances correspond to the S-matrix poles in the complex energy plane close to 
the physical region. This is in contrast to some other descriptions where resonance 
parameters are introduced using K-matrix poles or Breit-Wigner formulae with some 
ad-hoc background in each channel. Let us here remark that the Particle Data Group 
[p]] has misplaced the a meson parameters found in ||. They have been referred to 
under the name" fo (400 — 1200) Breit-Wigner mass or K-matrix pole parameters" on 
page 363 of |] instead of being referred under the appropriate title "/o(400 — 1200) 
T-matrix pole a/s". Similar misplacements were also made on pages 392 and 393 
for the /o(1370). Studying the analytical structure of multichannel amplitudes we 
shall show the important role played not only by poles but also by zeroes of the S- 
matrix. Knowledge of the poles and zeroes enables us to give in some cases a simple 
phenomenological approximation of the T-matrix. Our results on scalar-isoscalar 
resonances, masses, widths, branching ratios, coupling constants as well as phase 
shifts, inelasticities and cross sections in the three channels will be discussed for our 
different solutions and compared with the available data. 

In Section 2 we describe the two new solutions in addition to those previously 
presented in 0. Section 3 contains an analysis of positions of the S-matrix poles 
in the three channel amplitudes. Section 4 describes the influence of the S'-matrix 
poles and zeroes on phase shifts and inelasticities near the /o(1400) resonance. In 
Section 5 we comment on the limited applicability of Breit-Wigner approach in the 
multichannel meson scattering. Evaluation of branching ratios in two and three 
coupled channels is presented in Section 6 and coupling constants are discussed in 
Section 7. In Section 8 we present phenomenological parametrizations of multichan- 
nel amplitudes. Summary and conclusions are given in Section 9. In Appendix A 
the full formula for the Jost function of our three-channel model is supplied. Ap- 
proximated formulae for the pion-pion S'-matrix element, especially useful in vicinity 
of the /o(1400) resonance, are given in Appendix B. 

2 New and former solutions 

In we have briefly presented our three-channel model of meson-meson scat- 
tering. This model describes simultaneously nine reactions in a unitary way. In 
addition to the rcrc and KK channels an effective aa channel is introduced in order 
to describe important 4tt production and rescattering processes. Here we use ex- 
actly the same notation and definitions as in JFj. We furthermore give, in Appendix 
A, the full formula for the Jost function D(ki,k 2 ,k 3 ). Its knowledge is sufficient 
to construct all the S-matrix elements and all the physical quantities which will be 
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discussed in this paper. In |7| we have obtained four three-channel solutions called 
A, B, C and D based on \ 2 fits to the txtx 5-wave isoscalar phase shifts and inelas- 
ticities, and to the KK phase shifts. The parameters of the separable potentials, 
for all the different solutions, have been obtained by fitting the results of the fully 
coupled channel calculations to the experimental data. If the interchannel couplings 
are switched off then the KK pair remains unbound in all four fits. In solutions B, 
C and D the KK potential is repulsive, so a KK bound state cannot be formed. In 
the fit A the KK potential is attractive but its strength is too small to form a bound 
state. Parameters of the KK interactions, found in [7] , were different from the KK 
parameters obtained in || where a KK bound system could be formed in absence of 
the interchannel interactions. We should, however, remember that in a different 
mi phase-shift solution, obtained from the data analysis on an unpolarized target 



by the CERN- Munich group fTEfl, was used. 



In pj we have used several sets of the Tin phase shifts obtained in [p] from the 



data taken on a polarized target by the CERN-Cracow-Munich collaboration [17 
Now, one can ask if the data of || rule out any set of potential parameters leading 
to a KK bound state. We have therefore repeated the x 2 fits performed in for 
the "down-flat" tttt solution by adding an additional constraint, namely that the 
KK potential is sufficiently strong to obtain a bound state in the uncoupled case. 
As a result of these studies we have obtained two new solutions E and F which will 
be described below. Similar studies could be done for the "up-flat" solutions C and 
D which have less good \ 2 values as seen in Table 2 of [0. 



Table 1: Separable potential parameters for the solutions E and F. Values of (3 and 
m 3 are given in GeV. 



parameter 


Solution E 


Solution F 


An,i 


-.26349 x 10" 3 


-.13678 x 10" 3 


All, 2 


-.18316 


-.17845 


A 22 


-.60400 


-.52087 


A 33 


.17703 x 10 2 


-.73962 


Al2,l 


.28776 x 10~ 4 


.89713 x 10" 6 


Al2,2 


.036838 


.048444 


Al 3 ,l 


-.34811 x 10~ 3 


-.14046 x 10~ 3 


Al3,2 


.55929 


.095244 


A 23 


-1.5951 


-.024116 


Pl,l 


.16355 x 10 4 


.31518 x 10 4 


Pl,2 


1.1052 


1.0712 




1.4960 


2.1224 


A 


.092701 


1.4958 


m 3 


.675 


.680 
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Table 2: Comparison of four solutions fitted to the "down-flat" data of [8]. Second 
and third rows indicate presence or absence of bound states in KK and aa channels 
when interchannel couplings are switched off. In the remaining rows the different x 2 
values are specified. Numbers of experimental points are indicated in parenthesis. 
The x 2 values result from fitting with reduced rj errors as explained in [7]. 



Solution 


A 


B 


E 


F 


bound KK 


no 


no 


yes 


yes 


bound aa 


no 


yes 


no 


yes 


xl (65) 


63.0 


61.2 


64.8 


65.4 


XIk (21) 


15.9 


9.7 


17.7 


24.7 


Xl (30) 


13.2 


12.9 


12.6 


10.5 


Xlt (H6) 


92.1 


83.8 


95.1 


100.6 


Xf (30) 


36.7 


29.3 


33.4 


33.6 


Xlt (H6) 


115.6 


100.1 


115.9 


123.7 



The resulting potential parameters for the solutions E and F are presented in 
Table |. For the solut ion E we can notice a particularly strong repulsive coupling 
constant A33 in the aa channel. So the aa interaction for the solution E will not 
form any bound state when the interchannel couplings are set equal to zero. This 
is in contrast to the solution B where such a state exists. The aa bound state is 
also present in the solution F which has in addition a bound KK state. In Table 
|2| we give the binding properties when couplings between channels are switched off 
and the different x 2 values of the four solutions A, B, E and F. The x 2 values 
for the solution F are not as good as for the solutions A and B but they are still 
quite reasonable. The nir phase shifts of solution F increase by almost 180° at 
E pa 1350 MeV since at this energy a very narrow resonance (r pa 0.5 MeV) is 
created about 10 MeV below the aa threshold. This resonance makes also very 
narrow structures in the energy dependence of the inelasticity and phase shifts. 
They cannot be, however, uniquely confirmed by the existing data, so we tend to 
treat the solution F as an interesting but not well experimentally confirmed example 
of a phenomenological set of separable potential parameters. The x 2 numbers of the 
solution E are better than those of the solution F. They are comparable to those of 
the solution A, however, the solution B has the best x 2 ■ Therefore in the following 
mainly the three solutions A, B and E will be simultaneously discussed. We shall 
furthermore present new physical quantities like cross sections, branching ratios and 
coupling constants which have not been discussed in paper [[^. 

In Fig. 1 we compare, for energies above 1350 MeV [aa threshold), the inelastic- 
ities for the solution E in the three channels. Both tttt and KK inelasticities show 
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a minimum with a different depth near 1600 MeV and 1650 MeV, respectively. The 
aa inelasticity, however, has two minima at 1475 and 1675 MeV. Second minimum 
can be related to an additional scalar resonance at about 1700 MeV (/ (1710) ). 
Appearance of this resonance in addition to /q(1400) is a unique feature of the so- 
lution E. No such state exists in solutions A and B as it will be shown in detail in 
the next chapter. As we can see in Fig. 2 this state does not produce any strong 
increase of the inr phase shifts which is often expected in the vicinity of a resonance. 
This fact will be explained in the next section. 

In Fig. 3a,b we show inelasticities and phase shifts in the aa channel for three 
solutions A, B and E. Lack of resonances above 1600 MeV for the solutions A and B 
is responsible for a smooth energy dependence of the corresponding aa inelasticities 
above 1500 MeV. A smooth decrease is also visible in the aa phase shifts above 1450 
MeV. Influence of the /o(1400) resonance on the aa phase shifts can be seen for the 
solutions A and B as a small maximum around 1450 and 1400 MeV, respectively. In 
the solution E a monotonous decrease is disturbed by two structures around 1425 
and 1625 MeV related to the / (1400) and / (1710) states, respectively. The 
curve in Fig. 3b will be discussed in the next section. 

3 Positions of poles 

A knowledge of the positions of the S'-matrix poles in the complex energy plane 
is important for a physical interpretation of resonances. Some information about the 
positions of poles was already given in Jf| for four solutions: A and B corresponding 
to the "down-flat" data of || and C and D corresponding to the "up-flat" data. Now 
we shall discuss in more detail three sets of poles, all corresponding to the "down- 
flat" data, namely the solutions A, B and E which were described in the previous 
section. We think that the general structure of the solutions C and D is similar to 
those of sets A or B. In all the solutions S'-matrix poles appear on different sheets 
of the complex channel momenta k\, &2, k%. The sheets can be classified according 

to the signs of Imkx, Jmfe, Imk^. For instance, the notation h means that 

lmk\ < 0, Imki < and links > 0. Positions of the most important poles for 
the resonances / (500) , /o(980) and / (1400) were given in Table 3 of T7J. Origin 
of resonances can be studied by a gradual decrease of the interchannel coupling 
constants. In this way, starting from the case where all interchannel couplings are 
present, we arrive to the uncoupled case and obtain a trajectory linking the positions 
of a given pole from the fully coupled to the corresponding fully uncoupled case. 
The structure of poles for the various solutions is significantly different. Positions of 
different S'-matrix poles without and with couplings between channels are precised 
in Tables |3| to ^| for solutions A, B, E and F, respectively. In the two last columns 
the pole-sheet specification and their labels are given. 

Let us recall that a given pole in the uncoupled channel splits into four poles 
in the coupled channel case since coupling of a given channel to any other channel 
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Figure 1: Energy dependence of inelasticity in the tttt , KK and a a channels for 
the solution E 
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Figure 2: Energy dependence of irir phase shifts for the solution E 
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Figure 3: Energy dependence for the solutions A, B and E of a a inelasticities a) 
and era phase shifts b) . The curve denoted by (5 is the contribution of the double 
pole at k 3 = —iflz for the solution E. 
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doubles the number of related poles. This is a consequence of the analytic structure 
of the Jost function D(k\, k 2 , k^) describing three coupled channels. Furthermore, 
one should remember that 

D(ki, k 2 , h) = D*(-kl, -k*, -k* 3 ), (1) 

therefore to a given zero of D{k\, k 2 , k%) at k\, k 2 , k% there always exists a twin zero at 
— k\, — k%, —^3- In the complex energy plane a pole and its twin are symmetric with 
respect to the imaginary energy axis. In an uncoupled channel case zeroes and poles 
of S'-matrix lie symmetrically with respect to the real energy axis. This symmetry 
is broken when the interchannel couplings are switched on. Therefore asymmetry in 
the localization of the zeroes is crucial to understand the energy dependence of phase 
shifts and inelasticities in all the coupled channels. Now we shall discuss specific 
features of the different solutions. The positions of ^-matrix poles corresponding to 
the solutions A, B, E and F will be given in Tables 3, 4, 5 and 6, respectively. The 
"with couplings" column corresponds to the fully coupled model fitted to data while 
the "no coupling" column results from the same solution with interchannel coupling 
switched off. 

3a) Solution A 

For the solution A the tttt channel poles at E — (658 — z607) MeV and E = 
(1346 — z275) MeV in the uncoupled case evolve differently when the interchannel 
couplings are switched on (see Table |3|). The first pole leads to a set of poles related 
with / (500) at lower energy (poles I, II, and IV) while the second one splits into 
four states related with / (1400) at higher energy (poles VI to IX). The pole in 
the KK channel at E = (881 — z498) MeV, lying far from the physical axis in the 
uncoupled case, moves to the KK quasibound state at E = (988 — z31) MeV on 
sheet — h + (pole XIV). This is the /o(980) resonance which lies quite close to 
the physical axis and therefore it strongly influences behaviour of the tttt scattering 
phase shifts near the KK threshold. The pole trajectory linking the corresponding 
poles in the uncoupled and coupled cases is drawn as a solid line in Fig. 4a. In 
this figure we have also indicated at several intermediate positions the percentage 
strength of the reduced interchannel couplings. The same broad KK pole can also 
move to the pole XIII at E = (1038 - i204) MeV on sheet - + -. Its trajectory 
is drawn as a dotted line in Fig. 4a. Poles III, V, X and XI, lying on the real axis 
of energy below the tttt threshold, appear when the interchannel coupling strengths 
are large enough. Their influence on the ttit phase shifts is, however, small. 

The extremely wide pole in the third channel at E = (118 — z2227) MeV goes 
to very distant states at energies above 3000 MeV. These states, present only in 
solutions A and B, are quite model dependent. They cannot be visible in the energy 
dependence of phase shifts or inelasticities in the three open channels. We only 
give here their positions in order to illustrate the analytical structure of the model 
amplitudes. In many phenomenological applications such distant poles are treated 
as background described by additional parameters fitted separately to the data. We 
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Table 3: Positions of S'-matrix poles for the solution A (in MeV) 



channel 


no couplings 


with couplings 


Sign of 






ReE 


ImE 


ReE 


ImE 


Imk n Irnkx Imk a 


No 








564 


-279 




I 








518 


-261 


h + 


II 


TTTT 


658 


-607 


211 


o 


- + - 


III 








532 


-315 


h 


IV 








235 





+ + - 


V 








1405 


-74 




VI 


TTTT 


1346 


-275 


1445 


-116 


- + + 


VII 








1424 


-94 


1 


VIII 








1456 


-47 


h 


IX 








170 





H 


X 








159 







XI 


KK 


881 


-498 


418 


-10 


- - + 


XII 








1038 


-204 


- + - 


XIII 








988 


-31 


- + + 


XIV 








4741 


-4688 




XV 


aa 


118 


-2227 


3687 


-2875 


- + - 


XVI 








3626 


-3456 


+ 


XVII 








3533 


-579 


+ + - 


XVIII 



should stress that in our model there is no need to introduce any kind of artificial 
background. Large values of interchannel coupling constants (Table 1 in 0) are 
responsible for large pole shifts as seen in Table || and Fig. 4a. This applies also to 
other solutions considered below, especially to solution B (see Tables [I], f| to ^| and 
Fig. 4b). 

3b) Solution B 

For the solution B (Table |j) there are two wide poles lying far from the physical 
axis in the uncoupled tttt channel. The first pole at (733 — i583) MeV splits into 
four states I to IV. Pole II on sheet — h + and the corresponding zero lie close 
to the physical axis. Such a pole has an important influence on the tttt scattering 
amplitude below 1 GeV and can be related to the /o(500) resonance. The second pole 
evolves from (999-i323) MeV into a pair of poles at (1441 -zl25) and (1430 -zl49) 
MeV both related to the /o(1400) resonance. The two other shifted poles V and VIII 
remain close to 900 MeV. A very broad pole at (2744 — zl698) MeV in the uncoupled 
KK channel leads to four poles IX to XII lying far away from the physical region. 
Their influence on the scattering amplitudes is negligible. In Fig. 4b four interesting 
trajectories in the k 3 complex plane are drawn. The aa bound state at 1284 MeV 
(Imk 3 > 0) can either evolve to the / (980) resonance (pole XIII) at (992 — z34) 
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Im k% 
(MeV) 



500 
300 
100 
-100 
-300 
-500 



XIII 

b) ioo %■* 

.*■■ 

XIV^^"^ 

100 % 


E = 1284 MeV 
' % 


E = 1236 MeV 

o "A 


XVI 

.*' 

/'XV 

100 % 



-300 



-100 



100 

(MeV) 



300 



Figure 4: Pole trajectories in the complex momentum plane as a function of the 
percentage of the interchannel coupling strength: a) for the solution A in the KK 
channel, b) for the solution B in the era channel. Roman numbers and energies are 
taken from Table 3 for a) and from Table 4 for b). 
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Table 4: Positions of S'-matrix poles for the solution B (in MeV) 



channel 


no couplings 


with couplings 


Sign of 






ReE 


ImE 


ReE 


ImE 


Imk n Irnkx Imk a 


No 








332 


-114 




I 


7T7T 


733 


-583 


511 


-266 


- + + 


II 








512 


-266 


- + - 


III 








332 


-115 


h 


IV 








900 


-13 




V 


7T7T 


999 


-323 


1441 


-125 


- + + 


VI 








1430 


-149 


- + - 


VII 








942 


-25 


+ 


VIII 








3670 


-2263 




IX 


KK 


2744 


-1698 


3664 


-2240 


h 


X 








3102 


-940 


+ - - 


XI 








3104 


-904 


+ - + 


XII 








992 


-34 


- + + 


XIII 


oo 


1284 





1421 


-54 


+ 


XIV 








956 


-36 


- + - 


XV 


ao 


1236 





1411 


-85 




XVI 



MeV (dotted line) or to the / (1400) (pole XIV) at (1421 - z54) MeV (solid line). 
Similarly the so-called quasibound state at 1236 MeV (Imk 3 < 0) can either go to 
the pole XV at (956 - z36) MeV closely related to the / (980) (dotted line) or to 
the / (1400) resonance (pole XVI) at (1411 - z85) MeV (solid line). 

3c) Solution E 

In Table ^ the positions of poles for the solution E are shown. The first pole 
at (542 — i307) MeV evolves (as for the solution B) to four wide states. One of 
them (pole II) at (533 — «254) MeV on sheet — h + lies close to the physical region 
and can be related to /o(500) . When the interchannel couplings are switched on, 
the second pole at E — (1473 — zl50) MeV, coming from the tttt channel in the 

uncoupled channel case, creates four poles V to VIII. The pole V on sheet 

and a zero related to the pole VI on sheet — h + are responsible for an increase of 
the 7T7T phase shifts around 1400 MeV (see Fig. 2). Looking at inelasticities r/ n and 
r] K in Fig. 1 one can see small bumps around 1400 MeV caused by the pole V. 

The KK antibound state at 441 MeV transforms into two states IX and X lying 
far away from the physical axes in the KK and era complex momentum spaces. 
The KK bound state in the uncoupled case lies very close to the KK threshold at 
E = 990 MeV and can evolve, when interchannel couplings are switched on, into 
two narrow states (poles XI and XII) on sheets — I — and — h +• Pole XII on sheet 
— h + lies closer than pole XI to the physical region and therefore can be related 
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Table 5: Positions of S-matrix poles for the solution E (in MeV) 



channel 


no couplings 


with couplings 


Sign of 






ReE 


ImE 


ReE 


ImE 


Imk^ Imkx Imk a 


No 








600 


-355 




I 


7T7T 


542 


-307 


533 


-254 


- + + 


II 








533 


-246 


- + - 


III 








600 


-354 


+ 


IV 








1421 


-79 




V 


7T7T 


1473 


-150 


1441 


-106 


- + + 


VI 








1428 


-104 


- + - 


VII 








1466 


-38 


h 


VIII 








328 


-7 




IX 


KK 


441 





336 


-8 


+ 


X 








978 


-46 


- + - 


XI 


KK 


990 





990 


-34 


- + + 


XII 








1703 


-271 




XIII 


aa 


1565 


-112 


1648 


-67 


- + - 


XIV 








1673 


-77 


+ 


XV 








1624 


-175 


+ 


XVI 



to /o(980) . The last resonance at (1565 — zll2) MeV evolves into four states at 
about 1700 MeV. As it can be seen in Fig. 2 the titi phase shifts do not strongly 
increase around 1700 MeV. This fact is related to the positions of two singularities 
of the S'-matrix having the strongest influence on the titt scattering amplitude. One 

of them is the pole XIII on sheet and the second is a zero also lying on sheet 

related to the pole XVI on sheet H . The zero and the pole lie on the 

same sheet, so they partially cancel each other and their influence on the energy 
dependence of the titi amplitude is unusually small. The zero lies closer to both 
real and imaginary energy axes than the pole. Therefore in Fig. 2 a flat energy 
dependence of the tttt phase shifts between 1600 and 1750 MeV is seen as well as 
a smooth increase above 1750 MeV where the influence of the pole becomes larger 
than the action of the zero. 

Another particular case of influence of S'-matrix singularities on scattering am- 
plitudes can be seen in Fig. 3b. A strong decrease of the aa phase shifts for solution 
E, very well visible near threshold at 1350 MeV, can be understood if we take into 
account all the singularities of the S33 element of the S-matrix. For this solution 
the value of parameter /3 3 is very small (j3 3 = 92.7 MeV) so the double zero of S 33 at 
k 3 = —ifis is the closest singularity to the aa threshold. This zero of S33 is a result 
of the Jost function double pole at k% = —i(3 3 as it can be seen in the analytical 
expression for J 33 given in Eq. (A9) of 0. One can calculate the aa phase shifts 
keeping only this pole contribution to J33. The result is drawn in Fig. 3b as the 
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dotted line denoted by Comparison with the full calculation shows the domi- 
nance of this singularity over a large energy range. The poles related to the /o(1400) 
and /o(1710) only slightly disturb the energy dependence of the a a phase shifts at 
about 1425 and 1625 MeV. This example shows that one cannot extract a complete 
information about resonances from the energy dependence of phase shifts. One 
furthermore needs to know the analytical structure of the S'-matrix singularities. 

Let us here note that in the Particle Data Table there is a state called /j(1710) 
in a mass range similar to that of our / (1710), with a J = possibility not excluded. 
It has been recently shown that this J = assignment is more favorable than the 
J = 2 one @. 

3d) Solution F 

In the solution F there are both KK and oa bound states in the uncoupled 
channel case (see Tables |]and|]). In the fully coupled case, the a a quasibound state 
with a width about 0.5 MeV generates a very narrow 180° jump of the irn phase 
shifts at about 1350 MeV. The appearance of such a not well confirmed resonance 
cannot be excluded since the precision of the existing data is rather limited and the 
X 2 values for the solution F are still acceptable. Apart from this narrow state, the 
solution F is similar to other solutions so in the further analysis we shall not study 
in detail properties of this solution. 



Table 6: Positions of S'-matrix poles for the solution F (in MeV) 



channel 


no couplings 


with couplings 


Sign of 






ReE 


ImE 


ReE 


ImE 


Imk-K Imkx Imk a 


No 








691 


-511 




I 


7T7T 


554 


-377 


528 


-255 


- + + 


II 








673 


-410 


- + - 


III 








658 


-379 


+ 


IV 








1387 


-81 




V 


7T7T 


1407 


-80 


1428 


-93 


- + + 


VI 








1367 


-79 


- + - 


VII 








1447 


-77 


- - + 


VIII 








486 


-5 




IX 


KK 


497 





396 


-25 


- - + 


X 








967 


-31 


- + - 


XI 


KK 


968 





993 


-42 


- + + 


XII 








250 







XIII 


oa 


341 





270 





- + - 


XIV 








1349 


-0.2 


+ 


XV 


oo 


1353 





1349 


-0.3 


- + + 


XVI 
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Comparing different solutions we have seen that similar relatively narrow re- 
sonances can emerge from very different poles in the uncoupled channel cases de- 
pending on the set of the interaction parameters. Therefore one would need more 
data of a better precision to disentangle between phenomenologically good solutions 
A, B and E. In the next section we shall further discuss other properties of the 
meson-meson scattering amplitudes. 



4 Influence of S'-matrix poles and zeroes on phase 
shifts and inelasticities near /o(1400) resonance 

Properties of the scattering and production amplitudes in the energy region near 
1400 MeV can be understood provided that we know the positions of the S-matrix 
singularities, especially the poles and zeroes close to the physical region. In the 
previous chapter we have seen that the S'-matrix has many poles lying on different 
sheets. Now we would like to choose the most important poles and zeroes which 
influence the phase shifts and inelasticities and determine the resonance parameters. 

The S'-matrix elements S^- = 1,2,3) can be written in terms of the Jost 
function of different arguments, for example 

5ii = D(-k 1 ,k 2 ,k 3 ) 
D(k 1 ,k 2 ,k 3 ) 

Expressions for some other matrix elements can be found in fl~9| . All the S- matrix 
elements are inversely proportional to the Jost function D(k\, k 2 , k 3 ) which has a 
zero on sheet at the channel momenta kid (i — 1, 2, 3): 

D{k ld , k 2d , k 3d ) = 0. (3) 
The corresponding energy E d in the complex plane is given by 

E d = 2^Jk ld 2 + mi 2 = 2^Jk 2d 2 + m 2 2 = 2^k 3d 2 + m 3 2 , (4) 

where m.j denote meson masses. Energies E d on sheet corresponding to the 

/o(1400) resonance are collected in Table for our solutions A, B, E and F. The 
knowledge of this pole position of the S-matrix is, however, not sufficient to describe 
the pion-pion phase shifts and inelasticities even at the energy closest to the pole. We 
should also know the zero of the numerator D(—ki, k 2 , k 3 ) on sheet + + + close to the 
physical axis. This zero is in turn related to the zero of the denominator D(k±, k 2 , k 3 ) 
on sheet — h +. The corresponding energies on sheet — h + are also given in Table 
[7]. We can notice that both real and imaginary parts of energy are shifted on sheet 

— h + if we compare them with the corresponding parts of energy on sheet . 

The values of phase shifts and inelasticities depend on the zero of the numerator 
D(—ki,k 2 ,k 3 ) and on the zero of the denominator D(ki, k 2 , k 3 ). This means that by 
measurements of the irn phase shifts and inelasticities we cannot uniquely determine 
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one single value of the resonance energy and one value of the resonance width related 
to the imaginary part of energy Strong interchannel couplings are responsible for 
the energy shifts of zeroes found on different sheets. 



Table 7: Energy positions of S'-matrix poles related to /o(1400) 







Energy ( 


in MeV) 




Solution 


sheet 


sheet 


sheet 


sheet 






- + + 


- + - 


- - + 


A 


1405 - i74 


1445 - ill6 


1424 - i94 


1456 - i47 


B 


1411 -z85 


1441 - il25 


1430 - zl49 


1421 - i54 


E 


1421 - i79 


1441 - il06 


1428 - zl04 


1466 - i38 


F 


1387-^81 


1428 - i93 


1367 -i79 


1447 - i77 



In vicinity of the /o(1400) resonance one can, however, built up an approxima- 
tion to the S-matrix elements in three channels. This will furthermore allow us 
to understand the role played by the four poles whose energy positions in different 
sheets are given in Table [7| for our four solutions A, B, E and F. 

We make an expansion of the Jost function D(ki, k 2 , k 3 ) near its zero on sheet 
at (kid, hd, hd) 

D(ki, k 2 , k 3 ) « (ki - kid)di. (5) 

Let us denote by k in (i = 1, 2, 3) the zero position of the same Jost function on sheet 
— h +. Then in the first approximation 

D(-k 1 ,k 2 ,k 3 )^(ki-k* ln )ci. (6) 

Here d\ and Ci are complex constants. The quantity 1/di is the residuum of the 

1/D(ki, k 2 , k 3 ) pole on sheet and 1/ci is the residuum of the pole on sheet 

— h +• Here we have used the property of the Jost function expressed by The 

twin zero can sometimes be closer than the pole on sheet to the physical 

region and therefore it can strongly influence the S-matrix element. The hit S- 
matrix element can be then approximated by 

Sn = £±f /„ (7) 
ki - ku 

where fi = Ci/di and pi = —k\ n . Similarly the KK and aa S'-matrix elements can 
be written as 

s ^ = D(k u -k 2 ,k 3 ) _ k 2 + p 2 
D(ki,k 2 ,k 3 ) k 2 -k 2d 

D(k u k 2 ,-k 3 ) k 3 +p 3 
D(ki,k 2 ,k 3 ) k 3 -k 3d 
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where f 2 and f% are complex constants, p 2 is the kaon complex momentum on sheet 

— I — and j?3 is the a momentum on sheet h for which the Jost functions in 

numerators vanish. 

We have numerically checked, by comparison with the exact values, that the 
phase shifts and inelasticities for the solutions A, B and E, in the region of the 
/o(1400) resonance (effective mass range between 1350 and 1500 MeV), are qualita- 
tively well described by the formulae (0) to (||). For the solution B the agreement 
is even quantitative since the percentage error is only of the order of 10% or less. 
This error is obtained if the coupling constants /, are calculated using the first order 
derivatives of the Jost function. The errors can be reduced further if we modify the 
complex factors fi, f 2 and fs taking into account the second derivatives of the Jost 
function as explained in Appendix B. The solution E is somewhat particular since 
the region of / (1400) is also influenced by a wide resonance / (1710) (see Table 
|5|). The phase shifts in the a a channel are additionally affected by the second order 
zeroes at ^3 = ±i/?3 since the parameter (3% is small in that already discussed 

in Section 3c. 



5 Limited applicability of the Breit-Wigner ap- 
proach 

We should note here that zeroes of the Jost function D(k\, k 2 , ^3) on sheets 

, — h+, — I — and h are in general different, so the following inequalities 

between the corresponding complex momenta hold: p\ 7^ k^, P2 7^ k 2 d and ^3 7^ k^d- 

In particular, looking at the Table we can notice that the poles on sheet h are 

shifted towards higher energy in comparison with poles on sheet . Also the 

corresponding width related to the imaginary part of energy is considerably reduced 

on sheet h for the solutions A, B and E. This fact has important consequences 

on the energy dependence of phase shifts and inelasticity parameters which will be 
different from those obtained using the Breit-Wigner form. 

Let us now compare the approximations to the diagonal S-matrix elements (0) 
to (^) and the Breit-Wigner multichannel formula for the transition matrix elements 

t bw = ^ MY 

In this equation M is the mass and T the width of a resonance, q, Cj are real channel 
branching ratios and cj = Ti/T, Tj being the partial decay width. These formulae 
are valid if the resonance pole dominates the transition amplitude in the physical 
region and if background can be neglected. The corresponding S-matrix diagonal 
elements expressed by 

S? w = l + 2ihT* w (11) 
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are then written as 



BW s — [M 2 + l M(2T l - p] 

" s-{M 2 -iMT) ' 1 ' 

Remembering that s = E 2 = 4(k 2 + m 2 ) we can approximate (|T^) by 

where 

= — - m 2 - i— , (14) 

k tN 2 = ^-m 2 + tf(2T t -T). (15) 

While formally flTB|) looks very similar to (|7|), (|8|) and (^), it is very different from 
them since k^ and ki^ have to satisfy the inequalities: Rekny < RekiD, \ ImkiN |< 
| ImkiD |. These inequalities follow from an obvious inequality that each partial 
width Tj is smaller than the total width V. In fact each Breit-Wigner Sa matrix 
element depends only on three real parameters M, V and Tj while the matrix elements 
©' ® anc ^ (ED depend on three complex parameters: one giving the position of the 
zero in the denominator, the second giving the zero of the numerator and the third 
parameter being /j. 



and 



The Breit-Wigner formula fllCf ) is very often used to analyze experimental data 
in order to obtain a mass and a width of different resonant states []]]. It is widely 
believed that it should provide the same resonance parameters independently on the 
reaction channel in which the resonant signal is detected. Our analysis, however, 
puts some limits on a practical applicability of the Breit-Wigner approach. The point 
is that in each reaction channel not only the S-matrix pole plays a significant role 
but also the accompanying zero. The pole is common for all the reaction channels 
but the zero is different in each channel and it cannot be located simply by giving 
one number corresponding to a branching ratio or a channel coupling constant. This 
point will be discussed in detail in the two next sections. Here we can make the 
following remark. If one applies the Breit-Wigner formula to the data analysis in 
a particular channel then one can obtains distorted resonance parameters which 
are to some extend "averaged" over the pole and accompanying zero parameters. 
This effect might explain a fact that some resonant parameters "measured" in one 
coupled channel can be different from those obtained in another channels if the 
phenomenological model applied in the data analyses is essentially restricted to the 
Breit-Wigner formula. 
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6 Branching ratios 



Branching ratios are important parameters of hadronic resonances. Usually 
they are obtained from experimental data applying the multichannel Breit-Wigner 
formula (|To|) with some background parametrization in a limited range of effective 
masses. In our model we do not use any arbitrary background parametrization nor 
the Breit-Wigner parametrization. We fully exploit the knowledge of the analytical 
structure of the ^-matrix. Some S"-matrix poles close to the physical range can be 
related to the scalar meson resonances as already discussed in Sections 2 and 3. 
From the position of a given pole in the complex energy plane we can deduce the 
mass of the resonance and its total width. A determination of the partial decay 
width in presence of two and three open channels is, however, a more complicated 
issue which we are going to discuss in this section. 

6a) Definitions 

Let us recall that our model satisfies the unitarity condition for the S matrix: 
S + S = 1 and that the diagonal matrix elements are parametrized as 

S n = V3 e 2 ^, j = 1,2,3, (16) 

where rjj inelasticities and 8j the channel phase shifts. The nondiagonal elements 
are related to the non-diagonal reaction T-matrix elements by 

Tji = ^:-^=S jh j,l = l,2,3, j^l. (17) 
- z ■ 'kjki 



Expressions for Sji can be found in |TjJ. The diagonal T-matrix elements read 

^• = ^-(^•-1), 3 = 1,2,3. (18) 

T-matrix elements satisfy the unitarity equations. For example, in the first channel 
above the aa threshold 

ImTu = k x \T n \ 2 + k 2 \T 12 \ 2 + k 3 \T 13 \ 2 . (19) 

The total cross section in the first channels reads: 

87T 

a tot = /mTn . (20) 
h 



The elastic cross section is given by 



a 



ii 



S7i\T u \ 2 (21) 



and the transition cross sections from the channel 1 to the channels j = 2 or 3 are 
expressed by 

^■ = ^|7y| 2 . (22) 
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These cross sections satisfy the equation: 



a{f — a{ x + cr 12 + a 13 . (23) 

In Fig. 5 all the tot cross sections corresponding to the solution B are shown. At 
low energies we see a huge peak with a maximum near 600 MeV. It can be attributed 
to a very wide scalar a meson. Near 1 GeV one notices a very deep and narrow 
minimum which is related to the /o(980) meson. The fact that this resonance is 
seen as a dip and not, as in the most cases as, a maximum of the total cross section, 
is due to the special value of the elastic tctt phase shift which goes through 180° 
slightly below 1 GeV. If it happens below the KK threshold the total cross section 
is equal to 0. The next minimum of the total cross section is near 1500 MeV where 
another resonance /o(1500) appears. Up to the KK threshold the tttt scattering is 
elastic but above the KK threshold the tttt to KK transition cross section becomes 
visible. The transition to the aa channel is clearly visible above the third threshold 
at about 1350 MeV. 

Numbers of meson pairs produced in each channel are proportional to the cross 
sections defined above so for example in the channel we can define three branching 
ratios 

^ = -S^ J = 1,2,3. (24) 

These branching ratios obviously satisfy 

fen + b 12 + 613 = 1. (25) 

Cross sections similar to a\j and the corresponding branching ratios can be defined 
in channels 2 and 3 by changing index 1 into 2 and 3, respectively. In this way 
one can obtain nine branching ratios 6y = 1,2,3). All these quantities as well 
as the cross sections are functions of energy or the effective mass. Below the third 
threshold where energy is smaller than 2m 3 but greater than 2rri2 the branching 
ratio matrix 6^ reduces to 2x2 matrix containing only four nonzero matrix elements 
^11,^12,^21 an d 6 2 2- If the T- matrix is approximated by the Breit-Wigner formula 
(10) then the branching ratios are given simply by 

C = 4 (26) 
From (26) we infer that the corresponding branching ratio matrix contains iden- 
tical rows. 

6b) Discussion 

We shall discuss in detail the behaviour of branching ratios for the solution B, 
plotted in Fig. 6. Energy dependence of branching ratios for solutions A and E 
is similar to that for solution B. We do not expect that the energy behaviour of 
the branching ratios for the "up-flat" solutions (solutions C and D in @) will be 
qualitatively different from that of the "down-flat" solutions A and B. 
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E (MeV) 

Figure 5: Energy dependence of cross sections for the solution B 
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At first let us discuss the behaviour of the "elastic" branching ratio 6x1 as a 
function of energy in the pion-pion channel (see Fig. 6a). At the KK threshold 
611 = 1 but it decreases enormously steeply with increasing energy to about 0.6 
(solution A and E), 0.55 (solution B) or 0.5 (2-channel "down-flat" solution of J7|) 
already at the energy close to 1 GeV. This steep decrease is related to the opening of 
the KK channel. Next fen steadily increases to a maximum close to 0.9 for the above 
solutions at the energy reaching about 1.4 GeV. At this energy the third threshold 
is open for the solutions A, B and E. For E >1.4 GeV fen decreases very fast to a 
minimum close to 0.2 for three cases A, B and E and to for 2-channel "down-flat" 
solution. Then it rises again especially steeply for the latter solution and slower for 
the solutions A and B. This minimum is closely related to the presence of /o(1400). 

The fei2 branching ratio is equivalent to the transition probability of two pions 
into a pair of two kaons. Below the third threshold (E < 2m^) &12 = 1 — fen, so its 
behaviour is completely determined by fen. The fe i2 rises very steeply from the KK 
threshold attaining a maximum near 1 GeV, then decreases to a minimum value at 
about 1.4 GeV (see Fig. 6a). At this energy, however, the third threshold opens 
and the transition to the 00 channel becomes possible as stated by (|25|) . An effect 
of the interplay between three channels is that fei2 rises above 1.4 GeV attaining a 
maximum and then decreases rather slowly with energy for E >1.5 GeV. 

Above the third threshold in the tttt and KK channels a nonzero fraction of the 
channel total cross section originates from the transition to the 00 state. In the titt 
channel the branching ratio fei3 forms a regular maximum at E = 1.48 GeV, 1.46 
GeV and 1.48 GeV for the solutions A, B and E, respectively. The width of this 
peak is about 110 to 120 MeV (see Fig. 6a). The position and the width of this peak 
corresponds very closely to the parameters determined for the total width / (1500) 
by the Crystal Barrel Collaboration in different channels [|l|. 

In the second channel the fe 2 i element describes the ratio of the KK to titt 
transition cross section to the total KK cross section. At the KK threshold the 
annihilation cross section (X21 tends to infinity and the elastic cross section is finite, 
therefore the fe 2 i coefficient is equal to 1. With increasing energy it decreases to a 
quite deep minimum close to at about 1.4 GeV. At higher energy it rises slowly 
again for all three solutions. The transition element &23 starts from at the 00 
threshold, reaches rather quickly a maximum and then drops smoothly (see Fig. 
6b). 

Two branching ratios in the third channel, namely 631 and 632, decrease from the 
large values at the 00 threshold. Their behaviour is characteristic for the annihila- 
tion reactions. On contrary, the elastic branching ratio rises monotonically with 
energy as shown in Fig. 6c. Comparing Figs. 6a, 6b and 6c one can see that the 
Breit-Wigner formulae (^) are never satisfied above 1 GeV even in vicinity of the 
well defined resonances like /o(980) and /o(1400). 

6c) Average branching ratios near the / (980) resonance 

The / (980) resonance lies very close to the KK threshold and this fact has a 
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E (MeV) 

Figure 6: Energy dependence of branching ratios fry for the solution B. Dashed line 
corresponds to j = 1, dotted line to j — 2 and solid line to j — 3. In a) % — 1 (tttt 
channel), in b) % = 2 [KK channel) and in fig. c) i — 3 (era channel). 
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very important consequence on the experimental determination of the partial widths 
of this resonance for the decay into irn and KK channels. The nil channel is open 
below the KK threshold so for this channel we can determine the averaged tttt 
branching ratio b\\ choosing a range of the tttt effective mass ranges centered at the 
resonance mass M s with a maximum energy value equal to M max . This value should 
be larger than M s by more than the resonance width. This width is still not well 
defined experimentally. In |I[ the full width is 40 to 100 MeV. In our three channel 
analyses it was between 60 and 70 MeV and larger than these numbers by 10 to 
30 MeV for two channel fits. Therefore we can choose the value M max = 1100 MeV 
as the upper limit in the integration of the tttt or KK effective mass distributions. 
The average nir branching ratio is thus defined as follows: 

< b u >= — — — / b u (E)dE. (27) 

2(M max — M s ) J2M S -M max 

The average KK branching ratio over the same energy interval is 

< b 12 >= 1- < fen > . (28) 
Since b u = below the KK threshold 

< hi >= ^T7 ttt / h 2 (E)dE. (29) 

2{M max - M s ) J2m 2 

The mass values M s corresponding to the /o(980) resonance position on sheet — h + 
are 993, 989, 992 and 990 MeV for two channel down-flat model and three channel 
models A, B and E, respectively After integration over energy we have ob- 
tained averaged < bi 2 > values of 0.191, 0.156, 0.170 and 0.158 for the above four 
solutions corresponding to different sets of interaction parameters. Above numbers 
should not be, however, compared with the experimental value 21.9 ± 2.4% quoted 
as T(KK)/ r(7T7r) + T(KK) in the previous editions of the Reviev of Particle Pro- 
perties (see for example |2(J). The reason is that the numbers written under the 
title 'T 7r7r /(r 7r7r + r^x)" in 0,0] are the values of inelasticity coefficient rj av , defined 
below, and therefore they are not equal to the partial branching ratios of /o(980) . 
This difference is important since the branching ratios defined by (|24j) depend not 
only on the inelasticity coefficient r] but also on the nir and KK phase shifts. 

Experimentally the i] av value was obtained from the data on the 7r~p — > K^K^n 
reaction using the relation: 

2/rnN a(7r+7r- — >/o(980) — ► K%K%) . . 

where a u = 7c/(6k 2 ) is the unitary limit of the 7r + 7r~ — > KgKg scalar isoscalar 
cross section with k being the Kg momentum in the K S K S cm. system. Let us 
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remark that from the unitarity condition, below opening of the third channel, one 
has r) = rji = i]2- Thus one can define 



(1 rEmax \ 

p 5— / dE T) (E) (31) 

and use E max wl.l GeV as a good representation of the upper experimental limit of 
the KgKg effective energy (see [0] for a discussion of experimental uncertainties). 



Experimental values of i] av obtained in [21], |22|, |23[ are equal to 0.67±0.09, 0.78±0.03 
and 0.81±o;q4 respectively. Our calculations give /overvalues: 0.487, 0.496 and 0.494 
for the solutions A, B and E, respectively. We have also calculated rj av values 0.494 
and 0.761 corresponding to the two-channel "down-flat" solution of |7]] and to the 
two-channel set 1 result of H], respectively. Here we should stress that the last 
number is a result of the fit to the data for the reactions n~p — > K~K + n and 



n + n — > K~K + p of while the four previous numbers around 0.49 are based 
on inelasticities measured in the 7r~p^ — > 7T + 7T~n ||. The difference between above 
two-channel fits is therefore due to use of different experimental data sets. Problems 
with normalization of the tttc — > KK cross sections have already been discussed by 
Morgan and Pennington in (see in particular Fig. 4 therein) and Bugg, Sarantsev 



and Zou in [25 ]. The values of inelasticity 77 from ||, which were later used in J7J, 
correspond to the elastic titt — > tttt reaction. Although the errors are quite large the 
77 values found in Hj near the /o(980) resonance are visibly lower then those of P~5[] . 
This difference, however, could be partially explained by a possible contribution of 
channels other than tttt and KK . 

This difference between the data can also be seen in the calculation of < &12 >• If 
we use the parameters corresponding to the set 1 of || then < 612 >=0.055. This low 
value is related to rather low mass of / (980) equal to 973 MeV and to its width 29 
MeV obtained in |§. In the fits described in [^] we have obtained the / (980) masses 
close to 990 MeV and substantially larger values of the / (980) width. Authors of 
II and have also improperly used other definition of "branching ratio" in their 
Eq. (52) and in Eqs. (33) to (35), respectively. Those equations gave essentially the 
averaged values of 1 — rf which should not be directly compared to r\ av . 

One remark is in order here. If all the transition amplitudes are dominated by 
a single Breit-Wigner resonance as in (|I0| ) then the branching ratio 2x2 matrix has 
to satisfy the following relations: 621 = &11 an d 622 — ^12- These relations are not 
satisfied above the /q(980) resonance energy (see Fig. 6). It means among others 
that the transition amplitudes cannot be simply described in terms of a single Breit- 
Wigner formula. 

In Fig. 6a,b we see that the element 621 is in general larger than 612 especially 
near the KK threshold where 621 tends to 1 while 612 goes to zero. Nevertheless 
there is rather wide range of energies between 1.15 and 1.4 GeV where both ratios 
are quite close each other. 

6d) Average branching ratios in the range 1100 - 1420 MeV 
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In |27 the mass range 1100 to 1420 MeV has been chosen to compare the ratio 
of 7T7T to KK pairs produced by incoming pions on a polarized target at about 18 
GeV/c. The KK branching ratio was 

We have calculated the average values of bu in this energy region obtaining 0.189, 
0.175, 0.166 and 0.167 for the two-channel "down-flat" and three-channel solutions 
A, B and E, respectively. These values are higher than the experimental result of 
27] . One should, however, remember that the experimental errors of phase shifts 
and inelasticities are large (see 0) so the theoretical errors of 612 are also large. 
For completeness we give the value 0.099 corresponding to the data set 1 of 0. In 
order to get this value and the other four values written above we used the ordinary 
averaging procedure 



1 



'12 



M„ 



mm Mmin 



b 12 {E)dE 



(32) 



and not the definition (|29|) which is applicable only near the KK threshold. Here one 
should once again remind the confusion concerning the definitions of the branching 
ratios as discussed in subsection 6c, namely the value B av = 16 ± 1% calculated in 



should not have been compared with the experimental branching ratio 6.413,' 



70. 



6e) Average branching ratios near the /o(1400) resonance 

In presence of three open channels the branching ratio matrix b has 9 elements. 
As already discussed in 6b, behaviour of these elements in three different channels 
nn , KK and 00 is shown in Fig. 6 for the solution B. For the solutions A and 
E the general shape of curves is quite similar although numerical values of various 
branching ratios differ. In ]7J we have found that the width of the scalar resonance 
called here /o(1500) varies between 90 and 180 MeV depending on the solution. Now 
we can choose the energy interval 1350 MeV to 1500 MeV to present the averaged 
values of the branching ratios in the form of 3x3 matrix. The elements of the last 
line of this matrix are averaged over the energies larger than the third threshold 
energy equal to Im^,. For the solutions A, B and E we obtain 



636 





127 





237 


061 





844 





095 


300 





407 





293 






608 





163 





229 





054 





900 





045 





529 





286 





185 



and 






604 





133 





263 





063 





813 





124 





120 





156 





724 



27 



In the KK channel 622 dominates over 6 2 i and 623. Here the probabilities of the KK 
to 7T7T or 00 transitions are quite small. In the 00 channel, however, for the solutions 
A and B there are strong transitions from 00 to nn or KK transitions which actually 
are comparable to the transitions from the irn channel to era channels. 

In |28[] the branching ratios for the /o(1500) decay into five channels, rnr , 7777, 
7777', KK and An, are given as 29, 5, 1, 3 and 62 %, respectively. The two main 
disintegration channels are tttt and Air. In this model the An channel is represented by 
the effective aa channel and we also obtain large fractions for the averaged branching 
ratios 6n and 613. If we calculate the ratios 613/611 exactly at 1500 MeV then we 
obtain numbers 2.4, 1.2 and 2.3 for the solutions A, B and E, respectively. These 
numbers illustrate the importance of the An channel in agreement with experimental 
result of ||. In |9| the ratio r = (B[f -> KK]/B[f -> tttt]) h/fa = 0.24 ± 0.09 
is calculated (here by fo we mean / (1500)). If we define the ratio 612/611 then we 
obtain values 0.20, 0.27 and 0.22 for the solutions A, B and E, respectively. These 
values are close to r. From the partial decay widths of the / (1500) given in [ [2~5| 
one can calculate T K -j^/Y^ ~ 0.10 ±0.05 which is smaller than r but still consistent 
within the experimental errors. We know, however, that extraction of the branching 
ratios from experiment is a difficult task as it is, for example discussed in [[3(J. We 
should mention, however, that the average branching ratios depend quite sensitively 
on the energy bin chosen in the actual calculation as it can be seen in Fig. 6a. 
Furthermore we see that the branching ratio 6i 2 , corresponding to the nn ^KK 
transition, is very small around 1420 MeV, close to the position of our / (1400) 
resonance poles. This is in qualitative agreement with the small number for the 
KK branching ratio (3 %) given in |28| . 



7 Coupling Constants 



Coupling constants are useful physical quantities related to the residuum of the 
S'-matrix at its complex pole value 



s R = M — iMT, 



(33) 



where M is a resonance mass and V its total width. We define the coupling constants 
gi by the formula 

9i9i 



i^fs^ lim 
An s^s R 



Sij(s) 



(34) 



The product of the constants CiCj appearing in the Breit-Wigner formula (|T0"| ) is 
proportional to gigy. 



CiCj 



2MTJsr~ 



An 



(35) 
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In this way the diagonal Breit-Wigner coupling constant corresponding to the matrix 
element Su given by flT2| ) reads 

SiBW = 2MTi ^H ) (36) 

An k iD 

where is defined by (III). One can notice that the knowledge of gf/Air in the 
Breit-Wigner approach is equivalent to a determination of the partial width IY The 
two other independent quantities are M and V . 

The coupling constants corresponding to the approximations (0 to |) are: 

T- = &iV^fi(kd+Pi)- (37) 
An 



They are complex and depend on six parameters in contrast to (|36|) where only three 
parameters appear. 

Using (^4|) we have calculated the values of the coupling constants for three 
solutions A, B and E at two resonances / (980) and / (1400) . They are shown 
in Tables ||| and |], respectively. We can notice that the dispersion of the values 
| <7i | 2 / 47r corresponding to the nir channel is the smallest one since we have mostly 
fitted the tttt data. The worst situation is in the aa channel where there are no 
data available. In Table || we can notice that the KK coupling constants at the 
/o(1400) resonance are much smaller than the irn coupling constants contrary to the 
/o(980) resonance as seen in Table §. The knowledge of the coupling constants is 
not sufficient to describe fully the S'-matrix elements since the coupling constants 
are only related to the poles. We need also positions of the zeroes as discussed in 
Sect. 4. 



Table 8: Coupling constants of /o(980) at the pole — h + (in GeV 2 ) 



Solution 


l-nl 2 


Atx 


A 


0.37 


1.84 


B 


0.41 


1.33 


E 


0.40 


1.94 



8 Phenomenological parametrizations 
of multichannel amplitudes 

Following our studies described in previous chapters we shall here make some 
proposal for a simple parametrization of multichannel amplitudes. This parametriza- 
tion, based on existence of poles and zeroes of the 5*- matrix, has limited energy range 
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Table 9: Coupling constants of / (1400) at the pole (in GeV 2 ) 



Solution 




\92? 
47T 


4-7T 


A 


0.580 


0.106 


0.091 


B 


0.607 


0.233 


0.433 


E 


0.579 


0.107 


0.217 



of applicability. If there is some evidence of the presence of a well separated res- 
onance in some energy region one can use expressions ([?]) to (§J) to describe the 
diagonal S"-matrix elements. The corresponding diagonal T-matrix elements flT8] ) 
read: 

T« = f(/,-l + /,£^). (38) 

Ki Ki Kid 



From (jig) one can obtain the phase shifts and inelasticities in each channel 
provided that the three complex parameters (kid,Pi and fi), determined in a fit to 
some data, satisfy the unitarity conditions | Si | 2 < 1 in a given energy range. The 
off-diagonal S-matrix elements (|i~7|) can then be directly calculated as explained in 



19[] for the three channel case. 

The fitted parameters k^, related in all channels by the energy conservation 
(|J), provide us with the resonance energy position and its width. Other complex 
parameters Pi are, however, not related by (f|) since they correspond to zeroes on 
different sheets. If we consider a three-channel fit and if the resonance appears 

well above the third threshold then the pole lies on sheet which means that 

Imkid < 0. The effective parameters fi will take into account the second order 
derivative corrections to the approximated Jost function as explained in Appendix 
B (see in particular equations ( |B9"D to (|B12|) ). 

In Sect. 5 we have discussed the Breit-Wigner approximation. Our formula 
(j3~8"D can be reduced to the Breit-Wigner one if fi = 1 and Rek id ~ —Repi. Dy- 
namically fi = 1 means that the derivatives of the Jost functions in numerators 
and denominators of the Sa matrix elements (like in (0)) are identically the same 



see for example (|B10p). The phenomenology related to the formula (B8) is richer 



than that corresponding to the Breit-Wigner multichannel one ( |I~0| ) since it depends 



on six real parameters while the latter one depends only on three real parameters 
in each channel. If one adds the so-called elastic background to the Breit-Wigner 
amplitude Tjf then one gets the four-parameter formula 

T^^^ + e^T^, (39) 
where 5b is a background phase. Even introducing an inelasticity tjb to the back- 
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ground amplitude like 

Tii= VBG 7k~ 1 +r,Be2i5BT ^ w (4o) 

leads to a five-parameter formula. It is only when one enlarges the number of 
independent parameters in the T? w ([TOD , by allowing q to be complex, that (|40|) 
is essentially equivalent to our approximation (|3lf). The complex residues for broad 
overlapping states have been already introduced by several authors, in particular 
in |JT||. Limitations of the Breit-Wigner formula have been already pointed out 
in previous sections. Here we see that a knowledge of the S-matrix poles is not 
sufficient to construct the scattering amplitudes. We must also know positions of 
the nearby zeroes S'-matrix in order to describe the data with sufficient accuracy. 

A case in which one encounters more than a single resonance lying near the 
threshold should be treated differently especially if one wishes to study the energy 
range containing that threshold. Then in general more than one S'-matrix pole 
should be taken into account. In our analysis near 1400 MeV such a secondary pole 

appears on sheet h close to the physical region. Therefore this pole is included in 

Table [TO] showing average values of masses and widths of scalar resonances /o(500) , 
/o(980) and / (1400) obtained by us. Its mass is higher and its width is smaller 

than those of the pole. The three states /o(500), /o(980) and /o(1400) can be 

regarded as the model-independent resonances since they appear in all our solutions. 
Positions of other far lying poles above 3 GeV for the solutions A or B, of the pole 
around 1700 MeV in solution the E and of the pole at 1349 MeV in solution F are 
strongly model dependent. 

Table 10: Average masses and widths of resonances / (500) , / (980) and / (1400) 
found in our solutions A, B, E and F. Here errors represent the maximum departure 
from the average. 



resonance 


mass (MeV) 


width (MeV) 


sheet 


/o(500) or a 


523 ± 12 


518 ±14 


- + + 


/o(980) 


991 ±3 


71 ± 14 


- + + 




1406 ± 19 


160 ±12 




/ (1400) 


1447 ± 27 


108 ± 46 


h 



Ending this chapter let us remark that the formula ( p8|) can be used to approxi- 
mate the meson-meson transition amplitudes in different reactions where meson pairs 
appear. This would be the case in the meson production processes by hadrons, lep- 
tons and photons, in 77 reactions or in the central pp production, in heavy particle 
decays, in baryon-antibaryon annihilation and so on. 
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9 Summary 



We have presented two new solutions E and F resulting from fits to experi- 
mental data on ttti and KK phase shifts and inelasticities. Both new solutions are 
characterized by the existence of a KK bound state when the interchannel couplings 
are switched off. In the previously found solutions A and B the KK pair remained 
unbound in the uncoupled case (see Table §). In Tables |3| to |] we have given for each 
solution the set of S-matrix poles lying on different sheets. Knowledge of their posi- 
tions is necessary in order to evaluate masses, widths, branching ratios and coupling 
constants of scalar resonances studied here. We furthermore stress that not only 
poles of the S-matrix but also its zeroes play a very important role in evaluation of 
the resonance parameters. All zeroes are dynamically related. Also the knowledge of 
the pole trajectories as function of the interchannel coupling strengths is important 
to find origin of resonances. 

The fact that we have obtained different solutions while fitting the experimental 
data is mainly due to the limited number of data and due to their limited precision. 
We hope, however, that this situation can be improved in future when new data of 
better precision appear and the same model is applied to a wider class of reactions. 
Then one should be able to get nearly model independent information on the reso- 
nances lying sufficiently close to the physical energy axis. However, the poles lying 
far away from this axis will obviously be a subject of some model dependence. 

We have paid a special attention to the phenomenological analysis of the energy 
region near 1400 MeV where new scalar resonances appear [|] . A simple model of the 
S-matrix in three channels: irn , KK and 47r represented by an effective era has been 
constructed ((^) to @ and (|38|) ). This parametrization could be used to represent 
phenomenologically the meson-meson amplitudes in different reactions. The model 

is based on the knowledge of the pole on sheet common to all three channels 

and the zero which is specific for each channel. The diagonal .S-matrix elements 
depend on three complex parameters (positions of the pole and of the zero, and the 
strength). This is in contrast to the Breit-Wigner formalism where only three real 
parameters describe resonance amplitudes. The Breit-Wigner approach has limited 
applicability since it restricts too much the positions of zeroes and thus it does not 
take fully into account the interchannel dynamics. 

Branching ratios for different channels were defined and calculated for two- and 
three- channel models. We have studied the energy regions of the /o(980) and 
/o(1400) resonances and the region between them. We have noticed that different 
definitions of branching ratios were used in literature which lead to some confusion 
in past. The result is that there are no data on the mi and KK branching ratios 
of /o(980) and the incorrect determinations of them were no more given in the new 
edition of the Particle Data Tables 

Above the third threshold the branching ratios form a 3 x 3 matrix. Each row of 
this matrix describes three branching ratios in the particular channel, but only two 
of them are independent quantities (|25|) . Phenomenological information contained 
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in such matrix is much richer than in the Breit-Wigner approach where only three 
branching ratios can be used. We have studied energy interval 1350 MeV to 1500 
MeV and found small branching ratios corresponding to the transitions from the tttt 
to KK channel and vice versa. This is in qualitative agreement with findings of the 
Crystal Barrel Group ||| |9|. 



We have given formulae (^), flBTf ) and ( |3"7D for the coupling constants in the 
three-channel model and its different approximations. Calculated values of the cou- 
pling constants in the full model are reported in Tables ||and |[ The / (980) couples 
strongly to KK and the coupling of /o(1400) is much stronger for the tttt channel 
than for the KK channel. One should point out that the knowledge of the reso- 
nance mass and width and the coupling constants is not sufficient to give a good 
phenomenological description of the meson-meson dynamics. 

We have found that a fit of the tttt data of M together with KK phase shifts of 

[ p~5f| does require a /o(1406 ± 19) on sheet and a /o(1447±27) on sheet h 

(see Table 10). This seems to indicate that these data are quite compatible with 
Crystal-Barrel and other LEAR experiments which needs a broad /o(1370) and a 
narrower / (1500) Jl|. To make a more closer relation to the results obtained by these 
experiments one could perform a simultaneous analysis of the n~p^ — > 7r + 7r~n data 
and the data on pp —>■ 37r°, 57r° and other reactions. The model described by us can 
provide the production and scattering amplitudes of 9 reactions. These amplitudes 
can be applied in the description of the final state interactions between particles 
produced in different reactions. If needed our model of meson-meson interactions 
can be extended to describe more than three coupled channels (like f]7], rjr)' and so 
on). Our model is also sufficiently flexible to accommodate scalar mesons at energies 
higher than 1500 MeV, like the /j(1710) (see Table g). 

In the present version of our model the mass of the third channel is introduced 
as a free parameter. It is possible to construct a more involved model where one 
smears out this mass over a rather broad range according to the observed invariant 
mi mass distribution. Crossing symmetry constraints, if implemented, might also 
lead to reduce the number of parameters. Evidently new measurements should be 
performed near different meson-meson thresholds since even at the KK threshold 
the present data do not allow to obtain the /o(980) branching ratios. New data 
would help us to test different models of inter-meson interactions. It will lead us to 
a more profound knowledge of the structure of scalar mesons. 
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Appendix A 

Here we give the full formula for the Jost function D = D(k\, k 2 , k 3 ) of our 
unitary model built to describe the three channels: tctc , KK and the effective 2tt2tt 
(a a ). Using the same notation as in (0) we have 

D = d D + \ 00 I 01 D 1 - \ 02 I 22 D 2 - \ 03 I 33 D 3 , (Al) 

where 



d t 


— 1 Aji-Zjjj i — 0,1, 2, 3, 


(A2) 


Do 


= Dq — Ao2-^oi-^22C2 — ^03/01^33^3; 


(A3) 


c 2 


= A12 + /33(Al3A23 — A12A33), 


(A4) 


c 3 


— A13 + /22(Al2A23 — A13A22), 


(A5) 


D 


= d\d 2 d 3 — 2Ai2Ai3A23/n/22-^33 






— A 13 -/"ll -^33^2 — A^ 2 /n/22^3 — A 23/22-^33^1, 


(A6) 


D 1 


= — Ioid^Cl + A23(AnA23 — Ai 2 Ai3)/oi^22-^33 






— (Ao2-^00 + Ai2/oi)-^22C2 — ^O^OqI^C^^ 


(AT) 




= An + ^33(Ai 3 — A11A33), 


(A8) 



Di — Xoiil^Xn + — \n\jj)] + Ioo[d\dj — Xljluljj}} 

+ Aoj-^j[A23-^00^1 + AnA23-^oi + Ai2Ai3(Joo hi - / 2 i)] 

+Aij/oi«j + Xij\2zIoiIjji (A9) 



where i = 2, 3, j = 3 if i = 2 and j = 2 if i = 3. 

Appendix B 

We derive an approximated formula for the pion-pion S-matrix element given 
by (0). At first let us make an expansion of the Jost function in the denominator of 
Six near its zero on sheet at = (kid, k 2 d, k 3 d): 

D(k) = D(ki, k 2 , k 3 ) w (ki - hd)di + l -(h - k ld ) 2 si. (Bl) 

Here 

d _ dD(k 1 ,k 2 ,h) | k ld dD(ki,k 2 ,k 3 ) | ki d dD(ki,k 2 ,k 3 ) 
1 [ dki k 2d dk 2 k 3d dk 3 

and 



(B2) 
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Si 



d 2 D(k 1 , k 2 , k 3 } h/ d 2 D(k 1 ,k 2 ,k 3 ) k ld 2 d 2 D{k x , k 2 , k 3 ) 



+ 2- 



dk 2 
k ld d 2 D 



+ 2- 



k 2d 2 
k w d 2 D 



dk 2 



k 3d 2 



dk 2 



+ 2- 



<ld 



o 2 d 



k 2d dk 1 k 2 k 3d dhk 3 k 2d k 3d dk 2 k 3 
In ( |B~2| ) and ( [B~3| ) we have used relations 



(B3) 



k=k d 



k 



Id 



k 2 - k 2d « (ki - k ld )- — 
k 2d 



(B4) 



Id 



h - hd ~ (h - ki d ) — 
k 3d 



(B5) 



which follow from the energy conservation (Q). Similarly we make a power expansion 
of the Jost function in the numerator of S\i. 

Let us denote by k n = (ki n , k 2n , k 3n ) a zero of D(k\, k 2 , k 3 ) on sheet — h +: 



D(k ln ,k 2n ,k 



3n) 



0. 



(B6) 



Then using (ffl) we expand 



D(-kx, k 2 , k 3 ) « (fci - fci n *)ci + -{ki - k ln *) t 1: 



(B7) 



where 



ci 



dD(-k u k 2 ,k 3 ) k ln dD{-ki,k 2 ,k 3 



dk x 

k ln dD(-k u k 2 ,k 3 ) 



k 



2/i 



<% 2 



J k= 



and t\ is given by an equation similar to that for si in 



(B8) 



With the help of ( Bl|) and (|B7|) the pion-pion S'-matrix element reads 

, k\ — k\ n 



S 



ii 



ki — k 



id 



-./: 



i- 



where 



ft 



1 + 



"In 



1+ ill 



2 <i 



- A; 



id J 



(B9) 

(BIO) 
(Bll) 
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The function F(ki) can further be approximated by a constant if we limit ourselves 
to values k\ close to Rekid in the denominator of F(ki) and close to Rek\ n in its 
numerator: 

1 + ^ilmktn 

Similar expressions can be derived for the diagonal KK and era S-matrix ele- 
ments. 
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